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ABSTRACT 

Deep polarimetric Westerbork observations of the galaxy cluster Abell 2256 are presented, covering a frequency range of 325- 
377 MHz. The central halo source has a diameter of the order of 1.2 Mpc (18'), which is somewhat larger than at 1.4 GHz. With 
a = -1.61 ± 0.04, the halo spectrum between 1.4 GHz and 22.25 MHz is less steep than previously thought. The centre of the ultra 
steep spectrum source in the eastern part of the cluster exhibits a spectral break near 400 MHz. It is estimated to be at least 5 1 million 
years old, but possibly older than 125 million years. A final measurement requires observations in the 10-150 MHz range. It remains 
uncertain whether the source is a radio tail of Fabricant galaxy 122, situated in the northeastern tip of the source. Faraday rotation 
measure synthesis revealed no polarized flux at all in the cluster The polarization fraction of the brightest parts of the relic area is less 
than 1%. The RM-synthesis nevertheless revealed 9 polarized sources in the field enabling an accurate measurement of the Galactic 
Faraday rotation (-33 ± 2 rad m - in front of the relic). Based on its depolarization on longer wavelengths, the line-of-sight magnetic 
field in relic filament G is estimated to be between 0.02 and 2 /iG. A value of 0.2 yuG appears most reasonable given the currently 
available data. 
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1. Introduction 

Abell 2256 is one of the most massive clusters in the nearby 
universe. It has been studied extensively at X-ray, optical, and 
radio wavelengths. The galaxies associated with the cluster can 
be divided into three distinct populations, based on their kine- 
matics JBerrington et al.ll2002l) . X-ray images of the cluster re- 
veal sev eral substructures superimposed on a fai r ly smooth main 
source dSun et al.l l20()2l: iMolendi et all 120001: iRoettigeretal] 
119951; iBriel et al.lll991h . indicating that Abell 2256 is currently 
involved in a major merger event. 

The cluster has some interesting radio properties. It contains 
the largest n umber of head-tail galaxies of all known clusters 
JMilleret alJ b003). One of them has a straight, almost 1 Mpc 
long tail. There is a co mplex steep sp ectrum source in the east- 
ern p a rt of the cluster (|Rottgering et al. 1994; Masson & Maven 
Il978t lBridleetalJ |l979|). The cluster is embedded in a large 
diffu se, unpolarized radio halo JBridle & FomaloriBll976l ; iKimI 
1199 9; Clarke & EnBlin 200 ^ that dominate s the total flux at de- 
cametric wavelengths tCostain et al.l 119721) . The most striking 
feature is the large, bright relic source in the northwestern part of 
the cluster. Althotigh it is 20^0% linearly polarized at 1 .4 GHz 
(iBridle et al.ll979l;IClarke & EnBlinl2006l) . lJagersl(ll987h did not 
detect any linear polarization at 608.5 MHz, establishing a 3cr 
upper limit of 20% fractional polarization. 

There are strong indications that relic sources are caused by 
large scale structure formation (LSS) shocks compressing buoy- 
ant bubbles of magnetized plasma that have been emitted by ac- 
tive gala ctic nucle i in the past (Enfil i n & G o pal-Krishna 2001; 
EnBlin & Briiggenl 120021; iBruggenI 120031; iHoeft & BriiggenI 
2007,) . Shortly after their injection into the intergalactic medium 



(IGM), these bubbles have faded beyond detection limits due to 
synchrotron losses and adiabatic expansion. An encounter with 
a LSS shock wave would then compress the plasma, increas- 
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ing and aligning the magnetic field, and accelerate the electrons 
to energies enabling synchrotron emission at radio wavelengths. 
This model explains both the high fractional polarization that is 
usually observed in these relic sources and their peripheral posi- 
tion. 

De Bruyn & Brentjens (2005) have discovered several highly 
polarized structures in the direction of the Perseus cluster that 
were tentatively attributed to the cluster itself. Some of the struc- 
tures resembled buoyant bubbles in the IGM of the Perseus clus- 
ter, while others looked like long, straight shock fronts. One of 
these "fronts" was located at the western edge of the cluster at 
the interface to the Perseus-Pisces supercluster, exactly where 
such shocks are expected to occur. However, these structures 
have not been detected in Stokes /, hence the exact fractional 
polarization is unknown. The two most likely reasons for the 
non-detection of Stokes / are: 

- The sources are highly polarized and the Stokes / bright- 
ness is below the confusion limit of the Westerbork Synthesis 
Radio Telescope (WSRT) in these observations; 

- The structures have Stokes / structure at much larger scales 
than the Stokes Q and U structure, rendering it invisible at 
even the shortest interferometer baseline. 

The latte r situation is common in the Galactic synchrotron fore- 
ground ( IWieringa et al.l Il993h . A Galactic interpretatio n can 
therefore not be ruled out. However, iGovoni et alj (l2005h have 
observed similar structures in Stokes / near Abell 2255 with the 
VLA at 1.4 GHz. Abell 2256 was chosen to search for simi- 
lar objects because of its proximity, its dynamic nature, its relic 
source, and its moderate Galactic latitude of 32°. 

Wide field polarimetric images with noise levels well be- 
low 1 mJy beam"' are required in order to detect such sources. 
Additionally, the images can be used to find more polarized 
background sources in order to establish the Galactic contribu- 
tion to the Faraday rotation measure (RM) towards the cluster. 
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and to determine the decline of the fractional polarization with 
increasing wavelength, enabling a measurement of the thickness 
of the relic sources in RM-space. 

The spectrum of the di ffuse halo source has been rather un- 
certain. |Bridleietal] (|T973) favoured a spectral index (Sy <^ v") 
of a < -1.2 between 610 and 1415 MHz and a ~ -1.8 be- 
tween 151 and 610 MHz, based on maps where the halo was 
only marginally detected. More recently Clarke & EnBlin (2006) 
determined the 1 .4 GHz halo flux in more sensitive VLA D ob- 
servations, but could not derive the spectrum of the halo because 
of the large uncertainties in the halo flux measur emen t s avail able 
to them. They assumed that -1.7 < a < -1.25. lKiml(ll99 9') de- 
termined that a = -2.04+0.04, but this estimate was partly based 
on an erroneous value of the total cluster flux at 8 1 .5 MHz, and 
partly on maps in which the halo was only marginally detected. 
New measurements of the 35 1 MHz flux of the halo and the en- 
tire cluster are presented in Sect. |4] These values are combined 
with corrected literature values in order to obtain accurate spec- 
tra for the cluster as a whole and the diffii se halo source itself. 

The redshift of Abell 2256 is 0.058 (Berringtonet al.|!2002|). 
Assuming Ho = 71 km s"' Mpc"' (Spergel et al.. ,2003.) , 1' on 
the sky corresponds to 67.8 kpc at the cluster 

2. Observations 

The observations were conducted with the Westerbork Synthesis 
Radio Telescope, which consists of fourteen parallactic 25 m 
dishes on an east-west baseline and uses earth rotation to fully 
synthesize the uv-plane in 12 h. There are ten fixed dishes (RTO- 
RT9) and four movable telescopes (RTA-RTD). The distance be- 
tween two adjacent fixed telescopes is 144 m. 

The data set consists of two observing sessions of 12 h each. 
The first during the night of May 17/18 2004 with RT9-RTA = 
36 m and the second during the night of May 18/19 with RT9- 
RTA = 72 m. The distances between the movable dishes were 
kept constant (RTA-RTB = RTC-RTD = 72 m, RTB-RTC = 
1224 m). The uv-plane is therefore sampled at regular intervals 
of 36 m out to the longest baseline of 2736 m. Only the zero 
spacing is missing. The regular interval causes a grating ring 
with a radius of 80' at 350 MHz. At this frequency the -5 dB and 
- 1 dB points of the primary beam are at radii of 70' and 1 20' re- 
spectively. The observations are sensitive to angular scales up to 
90' at a resolution of 67" full width at half maximum (FWHM). 

The eight frequency bands were each 10 MHz wide and cen- 
tred at 319, 328, 337, 346, 355, 365, 374, and 383 MHz. The 
WSRT is equipped with linearly polarized feeds for this fre- 
quency range. The x dipole is oriented north-south and the y 
dipole east-west. The correlator produced 64 channels in all four 
cross correlations for each band with an integration time of 30 s. 
The on-line system applied a uniform lag-to-frequency taper A 
Hanning lag-to-frequency taper was applied off-line, effectively 
halving the frequency resolution. The pointing centre and phase 
centre were at the optical centre of the cluster at a = 17^03™45\ 
6 = +78°43'00" (J2000.0). 

The observations were bracketed by two pairs of calibrators, 
each consisting of one polarized and one unpolarized source. 
3C 295 and the eastern hot spot of DA 240 were observed be- 
fore Abefl 2256 and PSR B1937+21 and 3C 48 afterwards. 

The theoretical image noise in one Hanning tapered chan- 
nel of a single 12 h observation at 350 MHz is about 
1.56 mJy beam"' (uniform weighting). Band 1 was unusable on 
the second night due to a polarization calibration problem. Band 
4 suffered from a strange, broad band interference on the shorter 
baselines and was therefore discarded. Band 8 (383 MHz) was 



unusable in both nights due to man made interference. The four 
lowest and four highest channels of the remaining bands were 
discarded. Only the odd channels were selected because they are 
linearly dependent on the even channels due to the Hanning ta- 
per. In the remaining 28 channels per band, 60% of the data was 
usable. The expected thermal noise level in the integrated polar- 
ization maps is therefore 1.56 mJy beam"'/ V2 x 28 x 5 x 0.6 ~ 
0.12 mJy beam"' . This level can not be reached in the total inten- 
sity maps because they are confusion limited at appr oximately 
0.3 m Jy beam"' at this frequency and resolution (iMorgantil 
12004 . 

3. Data reduction 

Flagging, imaging, and self calibration we re performed with the 
AIPS-I-+ package dMcMulhn et al.l |2004|) . System temperature 
corrections, flux scale calibration, polarization calibration, iono- 
spheric Faraday rotation corrections, and deconvolution were 
performed with a calibration package written by the author and 
based on the AIPS-i-H- and CASA libraries. 



3.1. System temperatures 

The system temperature readings of the 36 m observing session 
were unfortunately unusable due to problems with the prepara- 
tion of the telescope for the observations. The system temper- 
ature readings of the 72 m session were usable, but had large 
spikes caused by RFI. The system temperature corrections were 
therefore only applied to the 72 m data, not to the 36 m data. The 
72 m system temperatures were filtered using a median window 
with a total width of 20 minutes, which was very effective in mit- 
igating the effect of short bursts of RFI. The difference between 
the X and y system temperatures of some bands in several anten- 
nas was of the order of 30-50 K, although for most band/antenna 
combinations it was less than 5 K. Because of the lack of system 
temperature data in the 36 m session, the 72 m data were used 
to determine the flux scale of the point source model used for 
self calibration. The 36 m data were later tied to the same flux 
scale in order to combine them with the 72 m data for the Stokes 
/ image. The 36 m data were not used for polarimetry because 
there was small, but significant polarization leakage even after 
cross-calibration. This may have been caused by temporal varia- 
tions in the difference in the system temperatures between the x 
and y receptors during the observation. 



3.2. Bandpass- and polarization calibration 

The flux scale and polarization leakages were calibrated simul- 
taneously using the unpolarized calib rator sources 3C 295 and 
3C 48. The Measurement Equation (ISault etal.lll996h for the 
visibility on the baseline i-j at a certain frequency is given by 



where 



V,, = 



< XX > < xy > 

< yx > < yy > 



(1) 



(2) 



represents the pure sky visibilities in terms of the cross correla- 
tions between the x and y receptors of the antennae in the base- 
line, J, and J; are the Jones matrices describing the properties 
of the antennae, and V,y is the observed coherency matrix. The t 
denotes Hermitian matrix transposition. 
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Table 1. Faraday rotation measures and position angles (north 
through east) at A^ - for the polarized calibrator sources. 



Source 



Rotation measure 
(rad m"- ) 



Position angle 

n 



PSR 1937+21 
DA 240 



+7.86 : 
+3.33 : 



0.20 
0.14 



291.4: 

122 : 



1.1 

3 



The fluxes of the calibrator source s were computed using the 
expressions and coefficients given i n iPerlev & TavloJ ([l_999u 
which extends thel Baars et al.l (Il977h flux scale to lower frequen- 
cies. Because each channel is individually tied to this flux scale, 
all sources appear in the images with their true spectral indices. 

The on-axis and off-axis polarization leakages at the WSRT 
are strongly frequency dependent with a period of the order of 
17 MHz. The leakage amplitude can increase from almost neg- 
ligible to 1.5% and decrease back to negligible over a frequency 
interval of approximately 8.5 MHz. It was therefore necessary 
to solve for all elements of the Jones matrices for each channel 
individually. The phases of the diagonal elements of the Jones 
matrix of RTl were fixed at 0. The remaining x-y phase dififer- 
ence. 



6n- = tan ' — + nn, 
U 



(3) 



where n is an integer, rotates Stokes U into Stokes V. It is for- 
tunately fairly constant across a 10 MHz band. Equation Q has 
only two unique solutions: one when n is even, the other when n 
is odd. One solution rotates the {U, V) vector to positive U, the 
other to negative U. In order to select the correct solution, one 
only requires the sign of the RM of the calibrator sources. The 
(sinusoidal) Stokes U spectrum is shifted by 90° towards smaller 
A^ with respect to the Stokes Q spectrum if the RM is positive 
and 90° in the opposite direction if it is negative. 

PSR 1937+21 and the eas tern hot spot of DA 240 both have 
a positive RM. According to iTsienI (|1982|) the RM of DA 240 
is +2.4 rad m"^ (no error quoted). According to A.G. de Bruyn, 
the RM of DA 240 plus a minimal ionosphere should be +3.7 + 
0.5 rad nT^. Using the low frequency front ends of the WSRT, he 
found that the RM of the pulsar, also without correcting for the 
ionospheric RM, is +8.5 ± 0.5 rad m~^ (private communication). 

The most accurate rotation measures for these objects to 
date are listed in Table [1] They were obtained after correcting 
for ionospheric Faraday rotation using the procedure outlined in 
Sect. 13.41 Because the shorter spacings were all alTected by solar 
interference, the polarization angle at each frequency was com- 
puted based on the average visibility of all spacings larger than 
700^. 

3.3. Phase calibration 

After the data were corrected for the more or less time indepen- 
dent effects described in the previous section, three phase-only 
self calibration iterations w ere performed on the 72 m data set 
JPearson & Readheadll 19841) . The sky model consisted of a list 
of 103 bright point source components with accurate positions 



' Note that the 327 MHz flux scale of WSRT observations has since 
1985 been based on a 325 MHz flux of 26.93 Jy for 3C 286 (the 



iBaars et alj dl977h value). On that flux scale, the 325 MHz flux of 
3C 295 is 64.5 Jy, which is almost 7% higher than the value assumed at 
the VLA and in this paper (A. G. de Bruyn, private communication). 



and /, Q, U, and V fluxes for each ch annel, supplem ented with 
a grid of 15525 CLEAN components (Hogbom'1974') in Stokes 
/. The CLEAN model was given a spectrum proportional to v ' . 
Both models were updated after every self calibration step. 
The initial models were obtained by: 

1. Making an image of a channel close to 350 MHz with all 
baselines > 600 A; 

2. Identifying all point sources brighter than 30 mJy beam"'; 

3. Solving for positions of all 103 sources in band 5 in the uv- 
plane (band 5 had very little RFI); 

4. Solving for /, Q, U, and V of all 103 point sources for each 
channel in the uv-plane; 

5. Solving for phases using the point source model and all base- 
fines >6QQA; 

6. Correcting the visibilities for these phases; 

7. Subtracting the point source model from the visibilities; 

8. Making (residual) Stokes / images of all usable channels and 
CLEAN them; 

9. Constructing the CLEAN model by averaging the CLEAN 
models of all channels. 

One self calibration iteration consisted of the following 

steps: 

1 . Adding the point source model to the visibilities 

2. Solving for phase corrections using the point source model 
and the CLEAN model on all baselines > 200 A; 

3. Solving for point source model fluxes; 

4. Subtracting the point source model from the visibilities; 

5. Making (residual) Stokes / images of all usable channels and 
CLEAN them; 

6. Constructing the CLEAN model by averaging the CLEAN 
models of all channels. 

The final CLEAN model was the average of the CLEAN 
models of individual channel maps that were produced after sub- 
tracting the point source component model. This lead to a very 
detailed and smooth representation of the extended emission in 
the CLEAN model. 

After the self calibration, data points were flagged if the 
absolute value of the residual visibility (corrected data minus 
model data) was larger than 5 times the root mean square (RMS) 
amplitude of the residuals for that band and baseline. 

3.4. lonosplieric Faraday rotation 

The ionospheric Faraday rotation at the WSRT is typically be- 
tween 0.2 and a few rad m"~. Daytime values of 5 rad m"^ are 
possible during solar maximum. A difference of 2 rad m"^ in 
ionospheric Faraday depth corresponds to a change in polariza- 
tion angle of > 100° at the lowest frequency (» 315 MHz). The 
polarized calibrator sources showed a difference of up to 20° in 
their polarization angle between the 36 m and 72 m sessions at 
the lowest frequencies, even though they were observed at the 
same time of the day on consecutive days. It was therefore nec- 
essary to compensate for the ionospheric RM in order to avoid 
de polarization. 

lErickson et al.l (l200 ll) showed that a simple ionospheric 
model fed with data from dual band GPS receivers installed at 
the VLA could predict the ionospheric rotation measure in most 
circumstances to better than 0.04 rad m"^. Sometimes, how- 
ever, they observed an unexplained discrepancy of the order of 
30° between the predicted and observed polarization angle of 
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Ionospheric Faraday rotation 
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Observed 



-4-202468 
Time since 2004/05/19 OOiOOh UTC [li] 

Fig. 1. Ionospheric RM as predicted by global GPS-TEC data 
and the US/UK World Magnetic Model. The dashed lines in- 
dicate the ±lcr predicted TEC/WMM rotation measures, in- 
terpolated between the 2 h intervals at which the data are 
published. The points with error bars are observed RMs of 
6C B165001.3 +791133, shifted up by 27.0 rad m'^. The hor- 
izontal axis is the time in hours since May 19 2004, 00:00 UTC. 



PSR 1932-1-109. This corresponds to a mismatch in RM of ap- 
proximately 0.6 rad nT^. 

Because there are only very few lines of sight through the 
ionosphere with the single dual band GPS receiver installed at 
the WSRT, the ionospheric Faraday rotation was computed using 
global GPS total ionospheric electron content (TEC) data and an 
analytical model of the geomagnetic field. The GPS-TEC data 
were provided by the Center for Orbit Determination in Europe 
(CODE) of the Astronomical Institute of the university of Bern, 
Switzerland. The geomagnetic field was computed using the 
US/U K World Magnetic Model (WMM) JMacmillan & OuinnI 
I2OOO1) . while the ionosphere was modelled as a spherical shell 
with a finite thickness and uniform density at an altitude of 
350 km above mean sea level. The geomagnetic field parallel 
to the line of sight was evaluated at the points where the line of 
sight from the WSRT towards Abell 2256 pierced through the 
model ionosphere. The TEC data are published for all even UTC 
hours. The expected uncertainty in RM towards Abell 2256, 
based on the TEC RMS uncertainty published by CODE, is typ- 
ically 0.07-0. 1 rad nT^ for each data point during the observa- 
tions. 

The ionospheric RM was also tracked using 
6C B165001.3 +791133, a polarized source at 53' from 
the phase centre. Its polarized flux is about 15 mJy at 350 MHz, 
yielding an apparent polarized flux of approximately 10.1 mJy 
after primary beam attenuation. The total Faraday depth as a 
function of time is 



<^tot(0 = l^src + 0ion(O, 



(4) 



where 0ion(O is the time dependent ionospheric Faraday rotation 
and 0src is the constant contribution of 6C B 165001. 3 +791133. 
Based on the thermal noise, 0tot could be estimated with an accu- 
racy of approximately 0.26 rad m"^ every two hours. This trans- 
lates to a Icr error in the polarization plane of approximately 13° 
at 324 MHz. 

Figure [1] compares the observed changes in ionospheric 
Faraday rotation with the predicted ionospheric RM based on 
GPS-TECAVMM data. The uncertainty in the GPS-TEC/WMM 



rotation measures is based on the RMS error in the TEC along 
the given line-of-sight. Uncertainties in the geomagnetic field 
are not incorporated. After discarding the outlier directly after 
sunrise, 0src = -27.0 ±0.1 rad nT^. 

The Abell 2256 visibilities were corrected for the iono- 
spheric Faraday rotation after instrumentally polarized point 
sources were subtracted in order to prevent residuals of these 
sources due to distortion of their PSF by the ionospheric Faraday 
rotation correction. 

3.5. Imaging 

The final dirty channel images in all Stokes parameters and the 
corresponding point spread functions (PSFs) were created from 
the corrected visibilities after subtracting the instrumentally po- 
larized point sources. The uv-plane was uniformly weighted. 
Because of a fractional bandwidth of 15%, the area under the 
main lobe of the PSF at the lowest frequency is 30% larger than 
at the highest frequency, hence a Gaussian taper was applied to 
the uv-plane before Fourier transforming in order to convolve 
the images to a common resolution of 67" FWHM. All maps 
are in North Celestial Pole (NCP) projection. The dirty maps are 
2048x2048 pixels large, with a pixel size of 15'.'6. The central 
1024x1024 pixels of the dirty ima ges were deconvolved using a 
Hogbom CLEAN (Hogbomlll974l) . 

While deconvolving the Stokes / channel images, CLEAN 
was constrained to only search for components at certain pixel 
positions (the mask). The mask was obtained by averaging all 
channel maps of a previous deconvolution, selecting all pixels 
where Stokes / was larger than some threshold, and expand- 
ing the mask to include all pixels where at least one of its 
eight neighbours was part of the mask. The threshold level was 
lowered after each self calibration iteration. The final thresh- 
old level for the mask was 1.5 mJy beam '. The loop gain 
was set to 0.3 and the CLEAN was stopped whenever either 
15000 iterations were performed or the maximum residual in 
the area where CLEAN components were allowed was less than 
0.5 mJy beam '. This is roughly one third of the thermal noise 
level in the individual channel maps. It was necessary to decon- 
volve so deep into the noise of individual channels in order to 
obtain the highest possible dynamic range in the average of 1 16 
channel maps. The CLEAN models were added to the residual 
images using a circular Gaussian restoring beam with a FWHM 
of 67". At this declination, the north-south to east- west ratio of 
the untapered WSRT PSF is 1.02. The use of a circular restor- 
ing beam is justified in this case because of the strong, circular 
uv-plane taper that was applied. 

After deconvolution, the channel maps were corrected for the 
total power primary beam of the WSRT, which can be approxi- 
mated by 



G(v, y) = cos kvy. 



(5) 



where v is the frequency in Hz, y is the angular distance from 
the pointing centre in radians, and k = 6.6 x 10"*^ s is a constant, 
namely the light crossing time across a 19.8 m aperture, which is 
roughly the effective diameter of a WSRT dish in this frequency 
band (25 m dish, 63% effective surface area). Finally, the best 
116 channel maps were averaged. 

The deconvolution of the Stokes Q and U maps was uncon- 
strained. The CLEAN was considered complete after 5000 iter- 
ations or if the largest residual on the inner 1024x1024 pixels 
was less than 1 mJy beam"'. These images were also corrected 
for the primary beam. Unfortunately the 36 m observing session 
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suffered from residual polarization leakage due to gain variations 
during the observations. It turned out to be difficult to obtain cor- 
rect gain solutions, and system temperature data were corrupted. 
Therefore the polarization images were created using only the 
72 m spacing. 



3.6. RM-synthesis 

The low frequency and high fract ional bandwidth of these ob - 
servations required RM-synthesis JBrentiens & de Bruvnll2005l) 
in order to avoid bandwidth depolarization. Images were created 
for Faraday depths of -800 to +800 rad nT^ in steps of 2 rad uT^. 

OfF-axis instrumental polarization was not corrected because 
it is limited to very specific Faraday depths after RM-synthesis. 
In the case of the WSRT the off-axis polarization consists of 
two components: a largely frequency independent offset and a 
frequency- and direction dependent oscillation with a period of 
17 MHz. After RM-synthesis, the offset ends up at (^ = rad m"- 
and the 17 MHz oscillation at \4>\ » 44 rad m"- (near 350 MHz). 
The sign of <p depends on the position angle with respect to the 
pointing centre. 

The noise level in the images ranges from 0.31 mJy beam"' 
near rad nT^ to the thermal noise level of 0.12 mJy beam ' 
for Faraday depths |0| > 300 rad m"^. The polarization images 
were corrected for the total intensity primary beam of the WSRT 
using Eq. (|5]l before the RM-synthesis was applied. 



4. Total intensity 

Figure |2] shows the central square degree of the final Stokes / 
image. The average frequency of the combined map is 351 MHz 
and the integrated bandwidth is 43 MHz spread over 52.5 MHz. 
The noise level near the centre of the map is approximately 
0.6 mJy beam"', dominated by residual calibration problems. 
The dynamic range of the full map (brightest sourcexentral 
noise level) is of the order of 3000:1. Sources are labelled 
according to Bridle et al.l (Il979l) . including the extension by 
iRottgering et all ( ll994i K 

Figure [3] displays a mildly super-resolved map of the same 
region. It was obtained by convolving the CLEAN model with a 
circular Gaussian beam of 26" FWHM. This m ap reproduces the 
filamentary structure seen in the relic area by (iClarke & EnBlinI 
[2006). The halo itself appears to exhibit filamentary structure 
too. 



Except for source AA, all sources mentioned in 
IRottgering et al.l (1 19941) were detected. For source AA, 
they list a total flux of 2.8+0.3 mJy at 1446 MHz. Although 
they mention the source has a flux of 8 mJy at 327 MHz, 
it is not visible in their 327 MHz map. H.J. A. Rottgering 
agrees that their flux estimate of source AA at 327 MHz is 
caused by an error in the flux estimation procedure (private 
communication). My non-detection gives a 3cr upper limit of 
1.8 mJy at 351 MHz. The source is also invisible in recent 
WSRT observations at 150 MHz, yielding a 3cr upper limit 
of 12 mJy (H.J. A. Rottgering, private communication). The 
specti-al index a (Sy oc y") between 351 MHz and 1446 MHz 
must therefore be larger than 0.3. 

The remainder of this paper focuses on the diffuse halo emis- 
sion that pervades the cluster, the steep spectrum source F, and 
the relic area containing filaments G, H, and the tail of source C. 
The complex involving sources A and B will not be discussed. 



4.1. Halo flux at 351 MHz 

The largest structure visible in Fig. |2] is the diffuse emission 
approximately centred on the X-ray source at J2000 position 
llH'""!' +78°37'55" (Ebelinget al. 1998), and extending to the 
sources F, O, L, H, and G. This halo source has been described 
bv several authors (iBridle & Fomalonli ll976t iBridle et al.|[l979l: 
iKimI 119991: IClarke & En6hn| |2006!)71t is considered to be flie 
main contributor to the total flux of the cluster at decametric 
wavelengths (Costain et al. 1972). 

The total flux of the halo at 351 MHz was estimated using 
the data selection from Fig.|2l The flux estimate is complicated 
by several bright, extended radio sources in the cluster, which 
have to be subtracted first. Because of their complex nature, the 
sources were subtracted in the image plane. Two different ap- 
proaches were used to estimate the halo flux in the area enclosed 
by the thick, grey contour in Fig.H] 

In the first approach the CLEAN model was set to zero in- 
side the areas enclosed by the white contours in Fig. |4] The 
model was subsequently convolved with a circular Gaussian of 
67" FWHM and added to the residual image. Let O be the set 
of no pixels enclosed by the grey contour and let fi be the set of 
«^ pixels enclosed by both the grey contour and a white contour 
The integrated halo flux is then the sum of all pixels in <i> that 



are not in fi, multiplied by 






and divided by the volume of 



the 67" circular Gaussian, which results in a flux of 696 ± 8 mJy. 
This is most likely a low estimate, because the sources A, B, C, 
D, F, and the eastern parts of G and H, are situated predomi- 
nantly in the brighter parts of the halo. The fainter south-eastern 
part of the halo therefore has a relatively large influence on the 
sum. 

In the second approach, the CLEAN model pixels in jj. were 
not set to zero, but were re placed by a radial basis function inter- 
polation d Carr et al]|200lh of the pixels that constitute the border 
of the areas enclosed by the white contours. This interpolated 
model was convolved with a circular Gaussian of 67" FWHM 
and added to the residual image. The result is shown in Fig. |4] 
The interpolation works extremely well for sources up to the size 
of source F The interpolation of the vast area of sources A, B, 
C, G, and H is less successful and appears a bit too faint in the 
area of A, B, and C, and a bit too bright in the area of G and H. 
The total flux was computed by adding the values of all pixels of 
the restored image that are in <I> and dividing by the volume of 
the 67" circular Gaussian. The result of 831 ± 10 mJy is prob- 
ably a somewhat high estimate because of the relatively high 
interpolated surface brightness of the relic area, which lacks the 
filamentary structure seen in the rest of the halo. The total halo 
flux is therefore estimated at 0.76 ± 0.07 Jy, which is the average 
of the two approaches. The uncertainty is mainly systematic. 

4.2. Radio spectrum 

The integrated radio flux of the entire cluster (halo, relic, and 
discrete sources combined) was determined from Fig.|2]by mea- 
suring the integrated flux within a circle centred at the opti- 
cal centi-e of Abell 2256 at J2000.0 position a = 17'•03'M5^ 
6 - +78°43'00". The radius of the circle was determined us- 
ing Fig. |5] The derivative of the integrated flux within a circle 
with respect to the radius of the circle is fairly high at small 
radii, where there is significant cluster emission. At radii > 5', 
the derivative drops sharply and settles at a more or less con- 
stant value at a radius of 10', which was therefore assumed to be 
the cluster radius. The integrated flux within 10' of the optical 
cluster centre is 3.51 + 0.06 Jy at a frequency of 351 MHz. 
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Fig. 2. Total intensity image of Abell 2256 observed with the WSRT in 2004. Contour levels are at (-4, -3, 3, 4, 5, 6, 7, 8, 9, 10, 
14, 20, 28, 40, 56, 80, 113, 160, 226, 320)x0.6 mJy beam^'. The beam is a circular Gaussian with a FWHM of 67". Crosses mark 
locations where point source components have been subtracted. The average frequency of this multi frequency synthesis map is 
351MHz. 



Table |2] lists flux estimates for the halo as well as the en- 
tire cluster including relics, head-tail galaxies, and background 
sources. All flu xes have been converted to the flux scale of 
iPerlev & Taylor! (11999) using the flux of Cas A, taking its secu- 
lar variation into account. A t frequencie s above 408 MHz, this 
flux scale is equivalent to the lBaars et al.l (11977) flux scale. 

No te that the flux at 8 1.5 MHz is twice the flu x from lBransonI 
(119671) . as suggested bv'Mas son & Maved (11978 *). The flux mea- 
surements from Kim ( 1999) have been read off Fig. 3 of that pa- 
per. The two measurements of the total cluster flux at 408 MHz 
and 1420 MHz appear to be systematically lower than all other 



cluster flux observations. It was therefore decided not to use 
these fluxes to determine the average spectra of Abell 2256 (halo 
plus relic and discrete sources) and its halo. 

The fluxes from Table |2] are plotted in Fig. |6] It is assumed 
that the cluster flux can be modelled as the sum of two spec- 
tral components, each with a constant spectral index between 
22.25 MHz and 2695 MHz: 



S(v)^S 



H.lOol 



\Q 



100 MHz ^ 



+ S 



R.lOol 



100 MHz/ 



(6) 
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Fig. 3. The average CLEAN model convolved with a circular Gaussian beam with a FWHM of 26". The contours are at (-4, -3, 
3,4, 5, 6, 7, 8, 9, 10, 14, 20, 28, 40, 56, 80, 113, 160, 226, 320)x91 juJy beam-i. 



where the first term is due to the halo and the second term is due 
to all other sources (relic and discrete sources). Fitting Eq. ^ 
to the total cluster fluxes, assuming uncertainties of 10% for all 
flux points with unknown errors, yields 



(7) 



This suggests that the spectrum of the halo may be less steep than 
the estimates by Co stain et al. ( 1972) ( -1.9 between 22.25 MHz 
and 8 1 .5 MHz) and lBridle et all d 19791) (- 1 .8 between 151 MHz 
and 610 MHz). 



5h,100 = 


6.0 ± 2.1 Jy 


an 


-1.65 + 0.22 


Srjoo = 


6.8 ± 1.8 Jy 


ffR 


-0.72 + 0.07. 



As Fig. |6] shows, the halo spectrum hes much closer to 
the ha lo fluxes derived in this paper and by 'Clarke & En61in| 
(|2006|) than to the flux estimates by Bridle et al. (1979) and 



iKimI (Il999l) . These sets of estimates are in fact mutually in- 
consiste nt. The reason is that the maps in this paper and the 
maps in lClarke & EnBlinI (l2006h ar e much deeper and dete ct the 
halo source over a larger area than lBridle et al.l ( 1197 9') and'Kim| 
(1999). The halo app ears to be ro ughly circular with a diameter 
of 12:2 at 1369 MHz dClarke & E nBlin 2006). The map in Fig.|4] 
shows that the halo is approximately square with a triangular 
south-eastern extension at 35 1 MHz. The sides of the square are 
18' long at 1 mJy beam"' and 13' at 3 mJy beam"'. The sides of 
the triangular extension are 6' long at 1 mJy beam"'. The halo 
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Fig. 4. The radio emission of the halo of Abell 2256 at 35 1 MHz. 
The contours are drawn at 1, 2, 3, 4, 5, and 6 mJy beam"' of 67" 
FWHM. In the areas enclosed by the white contours, the surface 
brightness was replaced by a radial basis function interpolation 
of the fluxes along the white contours. 



is only marginally detected in the 610 MHz map of iBridle et alj 
(1197 9) and at the noise level in their 1415 MHz map. They es- 
timated its diameter at 9' to 10'. Kim (199 9) estimates the halo 
size slightly larger ( 1 3 ' x 1 0'), but the halo is still only marginally 
above the noise in the 1420 MHz map presented in that paper It 
is therefore not surprising that the halo flux estimates in these 

papers were relatively low. 

Using the 1369 MHz halo flux of IClarke & EnBliiil ([2006h 
and the 351 MHz flux derived above, the spectral index of 
the halo is estimated at -1.5 + 0.2 between these frequencies. 
Including the 351 MHz and 1369 MHz halo fluxes in a joint fit 
with the total cluster flux points gives much better constrained 
values between 22.25 MHz and 2695 MHz: 



(8) 



The reduced x^ of the fit is 0.23, indicating that at least some 
fluxes in Table|2]are more accurate than stated. 



4.3. Source F 

Source F is a br ight, ultra steep spectrum source 
dMasson & MaveJ [1978; Bridle et al. 1979; Rottgerin Fet all 
|1994|) . An optical image from the Palomar Digital Sky Survey, 
overlaid with the Stokes / contours at 351 MHz is shown in 
Fig. [T] Source F has three components : Fl (southwest), F2 
(central), and F3 (northeast). iBridle et al.l (11979 ) suggested that 
because the spectral index between 1415 MHz and 610 MHz 
of the entire source is rather constant, Fl, F2, and F3 could 
be physically related. They suggested that the entire source 



'^H.lOO - 


5.9 ± 0.3 Jy 


an 


-1.61 ± 0.04 


Sr,\oo = 


6.7 ± 0.2 Jy 


ffR 


-0.72 ± 0.02. 




Fig. 5. The solid line is the total 351 MHz radio flux in Jy within 
a circle centred on the optical cluster centre as a function of the 
radius of the circle in arcminutes. The dotted line is the deriva- 
tive with respect to the radius. 

Table 2. Flux measurements for Abell 2256. 



V 

(MHz) 



(jy) 



C/H* Fif^ Reference" 



22.25 


100 ± 32 


C 


+ 


C72 


38 


41 ±6 


c 


+ 


W66 


81.5 


17 ±2 


c 


+ 


B67/M78 


151 


8.1 ±0.8 


c 


-1- 


M78 


351 


3.51 ±0.06 


c 


-1- 


B08 


408 


2.4 ± 0.24 


c 


- 


K99 


610 


2.165 ±... 


c 


+ 


B76 


1410 


1.13 ±... 


c 


+ 


B76 


1420 


0.75 ±0.11 


c 


- 


K99 


2695 


0.57 ±0.16 


c 


+ 


075 


2695 


0.666 ± . . . 


c 


+ 


H78 


151 


1.5 ±... 


H 


- 


B79 


351 


0.76 ±0.07 


H 


-1- 


B08 


610 


0.1 ±... 


H 


- 


B79 


1369 


0.103 ± 0.020 H 


+ 


C06 


1415 


0.02 ±... 


H 


- 


B79 


1420 


0.03 ± 0.009 H 


- 


K99 



" All values have been converted to the lPerlev & Tavlog i 19991) flux 
scale using Cas A. The error is omitted if no reliable error estimate is 
known. 

* Indicates whether the flux refers to the entire cluster (C) or the halo 
only (H). 

'' A "-" means that the point is for some reason not used to determine 
the spectrum of the halo and the cluster These reasons are described in 
the text. 

'' Re ference s: (C72) ICostain et al.' fi9lf); (W66) 'Willia ms et all 
( 119661) : (B67) lBransor] (11967): (M78) Masson & Mavei ( 1978); (B08) 
This paper; (K99) Kim ( 1999); (B76) Bridle & Fomalont ( 1976); (075) 
lOwen (1975); (H78) Haslam et al, U97&) : (B79) , Bridle et al., (1973): 
(C06) IClarke&Enfilinl JmoB) . 



is the radio tail of iFabricant et al.l (Il989h galaxy 122 at the 
northeastern tip of F3, and that F2 is a section of the tail oriented 
exactly along the line of sight. This wo u ld exp lain t he shell-like 
structure observed in 'Rottg ering et al.l (1 19941) and 'Mill er et aP 
(2003). Based on these maps, the size of F2 is estimated at 
45" X 90" a! 50 X 100 kpc. After subtracting Fig.gjfrom Fig.|2] 
in order to remove the contribution of the halo, the total flux of 
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Fig. 6. Spectra of the entire cluster and its two spectral com- 
ponents representing the halo and the rest of the cluster (see 
Eq. (|6]l). The labels indicate the source paper for the data points. 
Refer to Table|2]for details on the flux values. The solid lines rep- 
resent the sum of the two spectral components as computed by 
Eq. (|6]l. The dashed lines represent the individual spectral com- 
ponents. The long dashes represent the halo model and the short 
dashes the rest of the cluster. The thin curves show the spectral 
fit of Eq. (|6]l to the total cluster flux. The thick curves show the 
fit where the BOS and C06 fluxes of the halo have been used as 
extra condition equations on the spectral component of the halo. 
The total cluster fluxes of both solutions virtually overlap. 



source F at 351 MHz is 375 + 15 mJy. The 351 MHz flux of 
component F2 is 235 ±10 mJy. The peak brightness of F2 is 
185 mJy beam"' of 67" FWHM. 

Figure |8] shows the spectrum of component F2. The follow- 
ing data are included: 0.65 Jy at 151 MHz, 250 + 18 mJy at 
327 MHz, 235+10 mJy at 351 MHz, 97 + 6 mJy at 610 MHz, 
22 + 2 m Jy at 1415 MHz, a nd 23.6 + 2.5 mJy at 1446 MHz. 
Note that iBridle et al] (11979 !) do not mention an error for their 
151 MHz estimate. The spectrum appears to be considerably 
steeper at frequencies above 35 1 MHz than it is below that fre- 
quency, hence two power laws of the form 



S(v)^S 



100 



\ 100 MHz/ 



(9) 



were fit to the data. Between 151 MHz and 351 MHz 5ioo - 
1.06 + 0.07 and a = -1.20 + 0.05. Between 610 MHz and 
1446 MHz 5 100 = 2.15 + 0.40 and a = -1.71 ± 0.08. Based 
on the data in Fig. |9] it was determined that the spectral index 
of F2 between 338 MHz and 365 MHz is -1.2 + 0.2, which is 
consistent with the low frequency fit in Fig. [8] Extrapolating the 
35 1 MHz flux to 15 1 MHz using the 338/365 MHz spectral index 




Right Ascension (J3000) 



Fig. 7. Palomar DSS blue map of source F overlaid with 
351 MHz Stokes / contours. The contour levels are: 0.2, 0.28, 
0.4, 0.56, . . . mJy beam"', with a circular Gaussian beam of 26" 
FWHM. The role of Fabricant galaxy 122 is discussed in the 
text. 



yields 0.65 ±0.1 Jy, which is equal to the estimate bv lBridle et al.l 
(1979). 

The magnetic field i n F2 was estimated usin g the minimum 
energy formula given by lBeck & Krausd (l2005h . Unfortunately, 
the steepening in the spectrum compUcates the estimate. 
Approximating a convex spectrum with a single power law that 
is tangent to the spectrum at some point will overestimate the 
strength of the magnetic field because the synchrotron inten- 
sity, and therefore the cosmic ray energy density, is overesti- 
mated at both lower and higher energies than the frequency at 
which the power law is tangent. Because the flux is less aff'ected 
by synchrot ron losses at lower frequencies, the 151 MHz flux 
estimate by iBridle et all (1 19791) was used, yielding an average 
brightness of 0.17 ± 0.08 Jy arcmin"^ if a Gaussian shape of 
1 .7 ± 0.2' X 2 ± 0.2' FWHM is assumed. An injection spectrum 
of a X -0.7 was adopted. 

Based on its appearance at 1 .4 GHz (^ Rottgering et al.|[l994l : 
iMiller et al.ll2003l: IClarke & EnM in 2006) the most hkely mor- 
phologies are a spherical shell or a tube viewed along its axis. 
The line of sight through the source is approximately 100 kpc 
in case of a shell and a few hundred kpc in case of a tube. 
Given that the tail of a typical head-tail galaxy is somewhere 
between 100 kpc and 1 Mpc long, an estimated line of sight of 
500 kpc appears reasonable. The magnetic field in source F2 
was assumed to be tangled, whi ch seems valid given it s low 
polarization fraction at 1.4 GHz dClarke & EnBlinll2006l) . The 
minimum energy magnetic field in source F2 is estimated to be 
2.5 ± 0.3 X (Kq + l)i/(''+3) ^G for a Une of sight of 100 kpc and at 
most 1 .6 ± 0.2 X ( go + l)'^'°''^VG for a l ine of sight of 500 kpc. 

According to iBeck & Krausd (l2005h . a number density ra- 
tio of protons and electrons of Kq = 100 appears reasonable 
for several acceleration mechanisms, including Fermi shock ac- 
celeration, secondary electron acceleration, and plasma turbu- 
lence. Note that Kq is not the same as 7C, the energy density 
ratio of relativistic protons and electrons, which is traditionally 
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Fig. 8. Radio spectrum of source F2. The labels of the data 
points indicate its source (see Table |2]). The error bar on the 
151 MHz point indicates the uncertainty in the extrapolation 
of the 351 MHz flux with the 338/365 MHz spectral index 
to 151 MHz rath er than the uncertainty in the estimate by 
lBridleetalJ(fT979h . 



assumed to be 1 in a radio lobe plasma. Assuming Ko - 100 
gives 7.3 ±1.1 ijG and 4.8 + 0.8 yuG for Hnes of sight of 100 kpc 
and 500 kpc, respectively. These are the values used in the sub- 
sequent synchrotron age estimate. 

It is possible to estimate the age of a radio source based 
on the shape of its synchrotron spectrum. The shape of the 
spectrum reflects the history of the energy distribution of the 
injected relativistic electrons and the effect of various en- 
ergy loss mechanisms. In the following discussion it is as- 
sumed that synchrotron emission and inverse Compton scatter- 
ing of CMB photons dominate the energy loss, and that the 
energy distribution of the injected electrons follows a power 
law: tif^jE) dE cc E'^dE (e.g., Ginzburg & Syro vatskii 1965; 
Ivan der Laan & Perola|[T969l: iJaffe & Perolalll973h . It is further- 
more assumed that pitch angle scattering is effective, i.e., the 
momentum vectors of the relativistic electrons are distributed 
isotropically at all times. 

In the simplest scenario, the energy spectral index y and 
total power of the injected electrons is constant. In this case, 
the radio spectrum consists of a power law with spectral index 
ttiow — —{y ~ l)/2 below the break frequency and a power law 
with spectral index ahigh = -t/2 above the break frequency. 
The half life time of synchrotron emitting electrons emitting in 
a magnetic field with strength B in /vG is 



2.6 • 10'" years 



B2 



^rc 



B 



(l+z)vb 



(10) 




where Vb is the break frequency in MHz and Bic is the strength 
of a magnetic field with the same energy density as the CMB. 



17°06" 17°04°' 17°02"° 
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Fig. 9. Specti-al index map of Abell 2256 between 338 MHz and 
365 MHz. The Stokes / contours are drawn at 5, 10, 20, 40, 80, 
160, 320 cr of the 338 MHz map. The spectral index contours 
are drawn at -3 — hi in steps of 0.5. The accuracy is 0.46 at the 
40 sigma contour and 0.23 at the 80 sigma contour 



Assuming a CMB temperature of 2.725 K, its value is 3.238(1 + 
zf' juG. The power laws fitted to the data in Fig. [8] intersect 
at 400 + 180 MHz. Substituting this value for the break fre- 
quency and the minimum energy magnetic fields determined 
above yields an age of 51^^*^ million years if B = 7.3 yuG and 
77+27 million years if B = 4.8 /iG. 

The above model is somewhat problematic because it im- 
plies that y a: 3.4, whereas in most head-tail sources y » 2.4. 
It can therefore not be excluded that there is an additional break 
in the spectrum below 151 MHz, where the spectral index flat- 
tens to Q- X! -0.7. This situation can occur when the power of 
the injection of relativistic plasma has decreased sharply after 
a period in which it was high and steady. The resulting spec- 
trum has a spectral index of -0.7 below the hypothetical break 
at V < 151 MHz and a spectral index of -1.2 above this break, 
followed by an exponential decrease. This increases the age es- 
timates to > 83 million years if B = 7.3 jjG and > 125 million 
years B - 4.8 ^G. New observations below 150 MHz are needed 
in order to establish whether there is indeed an additional break 
in the spectrum. It is interesting to estimate the frequency range 
in which this break could occur 

Let us for the moment assume that F2 is indeed part of the 
tail of Fabricant galaxy 122. The galaxy is located at a distance 
of 2'6 (180kpc)fromF2in the plane of the sky. If F2 is indeed a 
section of the tail seen on-axis, one would expect that the galaxy 
travelled at least this far along the line of sight. The galaxy 
must therefore have travelled at least 0.25 Mpc. The current 
line-of-sight velocity of the galaxy is approximately 500 km s"' 
(iBerrington et al.ir2002l) . However, it is possible that the galaxy 
changed direction between F2 and F3 and is currently travelling 
approximately in the plane of the sky. Assuming that the actual 
velocity is of the order of the velocity dispersion of the cluster of 
approximately 1300 km s ' dBerrington et al.ll2002h . the elapsed 
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Fig. 10. The relic area after subtraction of the halo and the tail 
of source C. The lowest contour is drawn at 1 .5 mJy beam ' and 
all subsequent contours are scaled by V2. The resolution is 67" 
FWHM. 



time since starting the electron injection into source F2 is at least 
200 million years. A magnetic field of the order of 7.3 yuG then 
gives a break frequency of 26 MHz. This frequency will decrease 
if the source is older than 200 million years, but will increase if 
the magnetic field is lower than 7.3 /iG. Observations between 
151 MHz and 10 MHz will be possible with LOFAR in the very 
near future, but will be very difficult below 20 MHz. A complete 
spectrum from 10 MHz to 5 GHz will enable a more accurate 
minimum energy magnetic field estimate by properly integrating 
the cosmic ray energies over the observed spectrum, instead of 
using a single power law. Until a better constrained synchrotron 
age for this source is obtained, its association with galaxy 122 
remains uncertain. 



from the tail of source C to the border of the relic emission di- 
rectly south of source AC. There is a long, straight ridge south of 
the straight part of the tail of source C, at a position angle of ap- 
proximately 210° (N through E). In fact this source may extend 
north of the tail of C into the southern part of filament H. If that 
is indeed the case, the source is 8f8 (600 kpc) long. It appears to 
be unresolved at 26", limiting its width to less than 35 kpc. 

The flux of the relic area was determined after subtracting 
Fig. |4] from Fig. |2l Sources A, B, and C were excluded. The 
extended tail of source C was removed by interpolating fluxes 
of pixels nort h and south of th e tail using radial basis function 
interpolation (ICarr et all 1200 lb . The resuhing map is shown in 
Fig. [To] After correcting for the fact that the interpolated halo 
brightness in the reUc area in Fig. |4] appears to be on average 
1.5 + 0.5 mJy beam"' too high, the total flux in the relic area is 
1 .39 + 0.07 Jy. This is approximately a factor of three larger than 
the sum of the fluxes of the filaments G and H as determined by 
iRottgering et al.l (Il994l) . This discrepancy can be explained by 
the difference in uv coverage between the WSRT at 35 1 MHz 
and the VLA in B configuration at 327 MHz, which is not sen- 
sitive to s cales above 7^ Th e uv-coverage of the VLA-D obser- 
vation bv IClarke & EnBUnI (|2006) at 1369 MHz is much more 
similar to the WSRT at 351 MHz. Using their estimate of the 
relic flux, the spectral index between 1369 MHz and 351 MHz 
is -0.81 + 0.05. This is much st eeper than the 1446/327 MHz 
spectral index bv lRottgering et al.i ( 1994) , but sl ightly flatter than 
the ave rage 1703/1369 MHz spectral index that lClarke & EnBlinI 
(12006 ') found. 

The spectral index map in Fig.|9]shows that the 365/338 MHz 
spectral index is far from uniform across the relic area. The av- 
erage spectral index of filament H is -1.17 ±0.03, which is com- 
parable to its 1703/1369 MHz spectral index. The standard devi- 
ation of the spectral index distribution of filament H is 0.34. 

The spectral index of the brightest part of filament G is 
-0.70 ±0.1. Slightly north of that point is an east-west ridge 
with a 365/338 MHz spectral index of -0.3 to -0.4. There are 
several small areas with both steeper and flatter spectra. The 
average spectral index of filament G is -0.76 ± 0.03, which 
is significantly flatter t han its 1703/1369 MHz spectral index 
IClarke & EnBUnI (|2006|) . The standard deviation of the spectral 
index distribution is 0.37. 



4.4. Relic area 

The northwestern area of the cluster, dominated by the fila- 
ments G, H, and source C is perhaps the most intriguing part 
of Abell 2256. Most structure in that part is blended by the low 
resolution of Fig.|2] but Fig.|3]shows considerable detail. 

At the resolution of Fig. [3] the long, straight tail of source 
C is striking. It is 11^5 (780 kpc) long at 351 MHz, which 
is considerably longer than the 7.'08 (480 kpc) at 1446 MHz 
dRottgering et al.l 1 19941) . At approximately 6' from the host 
galaxy, the tail begins to bend: first towards the northwest, then 
just after source I to the west-southwest and finally back to the 
northwest near source S. The period of the tentative oscillation 
is approximately 6' (400 kpc) and the amplitude is about 0.'6 
(40 kpc). 

Although it is not as detailed as the 1369 MHz VLA C con- 
figuration image of lClarke & EnBlinI (12006). Fig.[3]clearly shows 
the filamentary nature of this part of the cluster The FWHM 
of the brighter parts of filaments G and H is of the order of 
of 60"±15" or 70±18 kpc. Filament G is approximately 10' 
(680 kpc) long and filament H is approximately 8' (550 kpc) 



5. Linear polarization 

Although the relic area is between 20% and 50% p olarized at 
1.4 GHz (IClarke & EnBhnll2006l; iBridle et al.|[T979l) . there was 
no linearly polarized emission at 35 1 MHz at rotation measures 
between -800 and +800 rad m"^ that could be attributed to 
the cluster. The relic area must therefore be substantia lly de- 
polariz ed. This is not surprising, as already noted by IJagersI 
(119871) at 610 MHz. A seai-ch in the area around Abefl 2256 
for structures similar to those a ttributed to the Perseus clus- 
ter (deBruvn & Brentiensll2005l) or found near Abell 2255 by 
iGovoni et al.l ( l2005l) did not uncover anything. The only visible 
polarized sources were the Galactic synchrotron foreground at (p 
between -20 and -24 rad m"^ and the collection of nine discrete 
sources (two doubles) listed in Table |3] 

These sources made it possible to estimate the Galactic con- 
tribution to the Faraday depth of the relic much more accurately 
than the current estimate of -4 ± 37 rad m"^ ( Clarke & EnBlinj 
2006). A coarse estimate is obtained by taking the mean value 
of the Faraday depths of the sources, which results in -29 ± 
3 rad m"^. 
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Table 3. Rotation measures and linearly polarized fluxes of discrete sources within 1°5 of the optical centre of Abell 2256. 



Nr.° 


Position* 







\p.L 


\Pt 




(J2000 a) 


(J2000 6) 


(rad m"-) 


(mJy beam ' ) 


(mJy beam ' ) 


1 


leHT'QTJ 


+79°06'16" 


-27.0 ±0.1 


10.1 


15.0 


2 


16H9"07;7 


+78°51'02" 


-28.0 ±0.8 


1.4 


1.8 


3 


16''51"45U 


+78°09'33" 


-29.9 ±0.8 


1.2 


1.7 


4 


17''00"3a;7 


+79°30'47" 


-20.8 ± 0.2 


5.5 


7.7 


5 


n'-oo'^ss'ie 


+79°30'32" 


-22.1 ±0.2 


6.1 


8.4 


6 


17''05"18';4 


+77°47'21" 


-24.7 ± 0.5 


2.0 


3.2 


7 


17''05"55^:5 


+77°06'30" 


-18.9 ±0.2 


6.5 


26.0 


8 


17h06"14M 


+77°07'47" 


-17.9 ±0.2 


7.1 


29.4 


9 


17''07"'17:0 


+78°06'31" 


-30.7 ± 0.5 


2.0 


2.5 


10 


17M8"44^:7 


+78°39'30" 


-49.1 ±0.3 


3.2 


4.2 


11 


17''23"'34^:7 


+78°29'18" 


-49.8 ± 0.5 


2.0 


3.4 



° Label of the source in Fig.llll 

* Peak of the polarized flux. 

'' Apparent polarized flux. 

'' Polarized flux corrected for primary beam attenuation. 



Galactic Faraday rotation 




n''zo' 



Right Ascension (J2000) 



le'so" 



Fig. 11. The Galactic Faraday rotation in rad nT^ overlaid on 
Stokes / contours. The Galactic RM is drawn at intervals of 
5 rad nT^ and the total intensity contours at -3, 3, 12, 48, 192 
xO.6 mJy beam '. The polarized sources from Table[3]are indi- 
cated with crosses. 



A better approach is to interpolate the Faraday depths of 
the sources. The result is displayed in Fig. [TT| The Galactic 
Faraday depth fiel d was estimated using a radial basis func- 
tion interpolation (ICarr et alj|200ll) . The interpolated Galactic 
foreground contribution in the direction of the relic is -33 ± 



2 r ad m ^. Subtracting this from the average RM of the relic 
that lClarke & EnBUnI (l2006h found (-44 rad m 2), yields a clus- 
ter contribution of -1 1 +2 rad m"^. 

Us ing an isothermal beta model determined by iMohr et al.l 
(11999 ') from ROS AT PSPC data (see Eq. (E)) and assuming that 
magnetic fields are frozen in (see Eq. (fT2Tl). this coiTesponds to 
a large scale magnetic field strength at the centre of the cluster 
of at least 0.08 jjG if the relic were situated at the same distance 
as the cluster centre, and at least 0.4 jiG if the relic is 500 kpc 
closer to the earth. 

The depolarization of the relic emission gives a handle 
on the conditions inside the filaments. Because there is no 
significant polarized flux in the relic area at 351 MHz, only 
upper limits to the polarization fraction can be established. 
Figure [12] shows a map of the upper limits to the polariza- 
tion fraction at 351 MHz. The maximum value of \F(ip)\ for 
- 8 00 ^ f ^ '^^P^ ^^'^ "^"^ ^^^ determined for each pixel. Refer 
to (Burn" 1966'; 'Gardner & Whiteoak|[T966t ISokoloff et ai][T998t 
Brentjens & de Bruvn 2005) for definitions of the Faraday dis- 
persion function F{(p). This map was then divided by the map 
in Fig. [2] The resulting upper limits are approximately a factor 
of four higher than upper limits obtained by dividing the image 
of the RMS value of \F{(p)\ in the same range of (p by Fig. |2] 
Because RMS - cr V/ if the noise distribution is Gaussian, this 
corresponds approximately to a 6cr upper limit. 

The fractional polarization is less than 1% in the bright- 
est par ts of filaments G and H at 351 MHz. The maps by 
Jageii (Il987|) give a 3cr upper lim it of 20% at these locations at 
608.5 MHz. lClarke & EnBlinl (1200 6) find about 28% polarization 
at both locations at 1369 MHz, with a resolution of 52" x 45". 
They furthermore find, without precisely specifying how it is 
spatially distributed, that the polarization fraction at 1703 MHz 
is generally about 8% higher than at 1369 MHz, which would 
imply roughly 30% polarization at these locations. 

Given the high degree of polarization at a resolution of 
the order of an arcminute at 1.4 GHz, the small spatial vari- 
ance in_^Mj_and_the smooth structure in polarization angle 
that IClarke & EnBlinl (|2006) find, it is unUkely that the sharp 
decrease in polarization fraction with increasing wavelength is 
due to beam depolarization. The naiTow channels, combined 
with RM-synthesis, rule out bandwidth depolarization at \<p\ < 
800 rad m -. This leaves differential Faraday rotation along the 



Brentjens: Abell 2256 at 350 MHz 



13 



351 MHz fractional polarization limits 



Depolarization of filament G 




IV'oe" 17''04°' 17'02" 

Right Ascension (J3000) 



iv'oo" 



Fig. 12. Upper limits (a; 6cr, refer to the text for details) to the 
polarization fraction in percent at 35 1 MHz at a resolution of 67" 
FWHM. Contours are drawn at 0.5%, 1%, 2%, 4%, 8%, 16%. 



line of sight, caused by the co-location of emitting plasma and 
Faraday rotating plasma as the most likely cause of the depolar- 
ization. The depolarization can be used to determine the Faraday 
thickness, or the extent in Faraday depth, of the radio emitting 
plasma. 

The depolarization data are plotted in Fig. [T3] Although an 
accurate determination of the Faraday thickness requires new, 
sensitive polarization observations between 0.4 and 1.4 GHz, 
it is interesting to obtain an estimate of the extent of the relic 
sources in Faraday depth using the available detections and lim- 
its. 

The uniform polarization angle structure at 1.4 GHz indi- 
cates that the magnetic field is relatively regular, perhaps even 
uniform at scales of a few arcminutes. If the relic has a uniform 
electron density, its structure in Faraday d epth can be approxi- 
mated by a uniform slab model (Burn 1966). The fractional po- 
larization as a function of A^ of a uniform slab is a sine function. 
The sine curve in Fig.[T3]corresponds to a slab with a full width 
of 21 rad m"^ in Faraday depth. As can be seen in Fig. [T3] it 
is impossible for a sine function to simultaneously satisfy the 
IClarke & EnBlinI (1200^) points as well as the upper limits pre- 
sented in this paper. The uniform slab model is therefore re- 
jected. The Gaussian that is plotted in Fig. [13] has a FWHM of 
4.7 rad nT^ in Faraday depth and satisfies all upper limits and the 
two points by Clarke & EnBIin (2006). A Gaussian was chosen 
because it is a function that smoothly goes from a high value to a 
low value and has an easy analytical Fourier transform. Because 
this particular Gaussian satisfies all points rather closely without 
exceeding the limits, the FWHM of 4.7 rad m"^ is considered a 
lower limit to the extent of filament G in Faraday depth. A more 
physical model will nevertheless have to wait until more sensi- 
tive observations between 400 MHz and 1.4 GHz are available. 
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0.001 
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■ "" ■ • .t Jageis 1987 



This work 
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2, 
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Fig. 13. Polarization fraction of the brightest part of filament G 
as a function of wavelength squared. The point bv'Ja gersI (Il987h 
is a 3cr upper limit. The point at the bottom right is half the limit 
from Fig.[T2]and is therefore equivalent to a 3cr upper limit. The 
horizontal bars indicate the A^ coverage of the observations. 

6. Magnetic field in filament G 

It is possible, with appropriate assumptions, to derive the mag- 
netic field strength in filament G from its depolarization proper- 
ties as displayed in Fig. [13] 

According to the isothermal beta model 



( Cavaliere & Fusco-Femianol 11976) determined by iMohr et all 
(ll999|) from ROSAT PSPC data. 






-3/J/2 



(11) 



where ne(^, is the thermal electron density at impact parame- 
ter b and line of sight distance between the cluster centre and the 

position / (see Fig.[T4]i, Wg.o - 3.60 + 0.06 x 10'^ h^^ cm"-^ is the 
central electron density, r^ = 342 + 3Qh:jl kpc is the core radius, 
andyS - 0.828 ± 0.06 is the ratio of the specific kinetic energies 
of galaxies and gas. The impact parameter of the brightest part 
of filament G with respect to the centre of the X-ray source is 6 '9 
(470±15 kpc). IClarke & EnBlinI (|2006|) established that the rehc 
must be on the front side of the cluster due to the small variance 
of the observed RMs. In the remainder of this section it is as- 
sumed that the brightest point of filament G lies approximately 
500 kpc in front of the centre of the cluster (I x -500 kpc). The 
expected thermal electron density in the neighbourhood of the 
filament is then approximately 4. 8 x 10" ^ cm"^. 



Following the reasoning of Mur gia et al.l (|2004|) . the mag- 
netic field strength is assumed to be a power law of the thermal 
electron density: 



IIBII 






(12) 



where a is the proportionality constant. If the magnetic energy 
density is proportional to the gas energy density, jj = /3/2. If the 
magnetic field is frozen in the plasma, then yu - 2y6/3. I only 
consider the latter case because the difference between the two 
is barely noticeable in view of the uncertainties involved. 

The relativistic electrons in relic sources such as filarnents G 
and H are assumed to be shock accelerated dEnBlin et al.lll998h . 
Merger shocks and accretion shocks typically compress the gas 
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Model electron densities in filament G 



Fig. 14. Coordinate system used for the computations in Sect. |6] 



by a factor of up to 4 JBlandford & Ostrikeilll978[) . For the sake 
of simpUcity it was assumed that the resuhing thermal electron 
density as a function of position along the line of sight has an 
offset Gaussian shape: 



■Je.thC/ - lo) = ne,{) 1 + (r - 1) e 



-4 log 2 



(13) 



where r is the shock compression ratio, w is the FWHM of the 
compressed area, and ne,o - 4.8 x 10""* cm"-' is the expected 
electron density near the relic as calculated using Eq. ( fTTI ). 

The synchrotron radiation is caused by relativistic electrons, 
which have a much lower density. The synchrotron intensity 
fror n a slice of plasma \y ith infinitesimal thickness dx is (see 
e.g. iPfrommer et al.ll2008l) 



«erell|B|P"" 

1(1) di oc ..':'':"' "_, d/. 



IIBIP + B- 



(14) 



ic 



where Bic = 3.238(1 + z)^ //G is the contribution from inverse 
Compton losses due to cosmic microwave background photons 
and ||B|| is given by Eq. (fT2] i. 

It is assumed that there are only significant amounts of rel- 
ativistic electrons in the compressed area. More specifically, the 
relativistic electron density is assumed to be a Gaussian with 
unit peak and the same FWHM w as the compression peak in the 
thermal electron distribution. The absolute number of relativistic 
electrons is unimportant for depolarization arguments. Only the 
shape of the distribution matters. Figure [15] shows an example 
of a thermal electron density profile with its accompanying rel- 
ativistic electron density profile. Note that the vertical scale of 
the relativistic profile is arbitrary. More precise modelling is not 
useful at this point because there are only upper limits available 
below 1369 MHz. 

The Faraday depth (f)(1) of a given location is obtained by 
integrating 



d(h 1 

^=0.81«x«J., 



(15) 



from / = -2w to the location of interest. The Faraday dispersion 
function 
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Fig. 15. Model electron densities used for estimating the mag- 
netic field in filament G. The model densities are plotted as a 
function of position along the line of sight with respect to the 
centre of filament G. The FWHM of the peaks is 100 kpc. The 
vertical scale applies to the thermal electron density only. The 
peak gas compression ratio in this plot is three. 



From symmetry arguments one expects the FWHM thick- 
ness of the emitting filaments in the relic to be of the order of 
30-100 kpc. Figure[T6]shows the peak magnetic field in the relic 
as a function of the shock compression ratio for FWHM thick- 
nesses from 25 to 400 kpc. The peak field strength was derived 
by finding the value of a that created a FWHM of F((p) of ei- 
ther 21 or 4.7 rad m"^, which was subsequently substituted into 
Eq. dni). 

As Fig. [T6]shows, the field estimate depends on the assump- 
tions about the shock compression ratio and the extent of the 
filament in both geometric space as well as Faraday space. A 
value of the order of 0.2 yuG appears reasonable assuming a relic 
thickness of the order of 30 kpc and an extent in Faraday depth of 
4.7 rad m"^ FWHM. However, due to the uncertainties involved 
this value may be off' by a factor of up to 10. 



7. Conclusions and future work 

The spectral index of the diffiise central halo of Abell 2256 
is -1.61 ± 0.04. Although t his value is steep er than the 
610/1415 MHz value of -1.2 in Brid le et all (IT979I) . it is some- 
what flatter than their 150/610 MHz spectral index (-1.8) and 
considerably flatter than the value of -2.04 derived by iKimI 
d 19991) . The reason for this is twofold: 



F(<p) = t(4>)/^- 



(16) 



- in IClarke & EnBlinI (1200 6') and in this work, the halo was 
de tected over a larger ext ent than in p revious observations 
bv lBridleetai] d 19791) and lKimI d 19991) due to the increased 
sensit ivity o f the observations; 

- Kim (1 19991) did not correct the 81.5 MHz fl ux JBransonl 
119671) as proposed bv lMasson & Maved (Il978h . 

The observed spectral index puts the classical minimum en- 
ergy magnetic field in the range of 2-4 fiG, and the hadronic 
minimum energy field between 4 and 9 juG JClarke & EnBlinI 
2006). 

Despite the deep polarization images obtained after RM- 
synthesis, no polarized emission was detected that could be at- 
tributed to Abell 2256. The complete depolarization of filaments 
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Fig. 16. Estimates of the magnetic field parallel to the line of 
sight in filament G. 



- apart from the relic area, there are simply no other bubbles 
or large scale shocks in Abell 2256 that emit radio waves. 

The first point can be tested by applying RM-synthesis to 
deep polarimetry at 1.4 GHz. Discriminating between the other 
two possibilities would be difficult at best in a single cluster 
One way to move forward is by observing a larger sample of 
nearby clusters at high Galactic latitude using RM-synthesis 
around 350 MHz. A fractional bandwidth of 20% is then suf- 
ficient to separate cluster emission from Galactic foreground 
emission. If RM-synthesis does not reveal any bubble-like polar- 
ized structures, similar to the ones towards the Perseus-cluster 
dde Bruvn & BrentiensI l2005h . it must be concluded that the 
structures detected towards the Perseus cluster really belong to 
the Galactic ISM. At this moment, however, there is not enough 
evidence to draw a final conclusion. 
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G and H was expected because of the relatively high density en- 
vironment in which they reside. The upper limit to the fractional 
polarization of less than 1 % is consistent with the 608.5 MHz 
upper limit by iJagersI (Il987l) and the depolarization between 
1705 MHz and 1395 MHz that'Clarke.&En61in| (E^ found. 

Reliable measu rements of the magn etic field in relic sources 
are relatively rare. iBagchi et al.l (Il998l) used inverse Compton 
scattering to derive the cluster scale magnetic field in the rehc 
of Abe ll 85. They found a value of 0.95 ± 0.10 ^G. IChen et all 
(|2008|) derived classical equipartition field strengths of 0.63 and 

1.3 juG for the relics in 0917-1-75 and 1401-33 respectively. 
Their 3cr inverse Compton lower limits are 0.81 and 2.2 fiG. 
lEnBUn et all (1 1998b listed an equipartition field of 3h'^['' juG for 
the entire relic in Abell 2256. Although the estimate of the line- 
of-sight field in filament G presented in this work is not very 
accurate (0.02-2 juG with 0.2 fiG in case of reasonable assump- 
tions about the shock compression ratio and path length along 
the line of sight), it is consistent with the quoted literature values 
for magnetic fields in relic sources in general. An actual mea- 
surement of the depolarization of filament G can improve the 
estimate significantly, although the major uncertainty is in the 
physical thickness of the relic along the line of sight, which is 
very difficult to obtain. 

Polarization observations covering frequencies between 

1.4 GHz and 400 MHz, where the fractional polarization de- 
creases rapidly, would make it possible to recover the (linearly 
polarized) emissivity of filaments G and H as a function of 
Faraday depth along the line of sight. The high fractional po- 
larization at 1 .4 GHz ensures that the magnetic field has no re- 
versals along the line of sight. Therefore, the Faraday depth of a 
certain location in the relic sources should be a monotonic func- 
tion of the geometric distance between the front of the sources 
and that location, enabling the first 3D reconstruction of an ex- 
tragalactic synchrotron source. 

The fact that no evidence was found for polarized emission 
from LSS shocks or buoyant bubbles further away from the clus- 
ter centre than filaments G and H could have several reasons: 

- they are depolarized due to internal Faraday dispersion, just 
like sources G and H; 

- there are no shocks oriented edge-on, only face-on, resulting 
in an undetectably low surface brightness; 
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